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Summary 

The properties of ceramic matrix composites (CMC’s) are 
known to display a considerable amount of scatter due to 
variations in fiber/matrix properties, interphase properties, 
interphase bonding, amount of matrix voids, and many 
geometry- or fabrication-related parameters, such as ply 
thickness and ply orientation. This paper summarizes prelimi- 
nary studies in which formal probabilistic descriptions of the 
material-behavior- and fabrication-related parameters were 
incorporated into micromechanics and macromechanics for 
CMC’s. In this process two existing methodologies, namely 
CMC micromechanics and macromechanics analysis and a fast 
probability integration (FPI) technique are synergistically 
coupled to obtain the probabilistic composite behavior or 
response. Preliminary results in the form of cumulative prob- 
ability distributions and information on the probability sensi- 
tivities of the response to primitive variables for a unidirectional 
silicon carbide/reaction-bonded silicon nitride (SiC/RBSN) 
CMC are presented. The cumulative distribution functions are 
computed for composite moduli, thermal expansion coeffi- 
cients, thermal conductivities, and longitudinal tensile strength 
at room temperature. The variations in the constituent proper- 
ties that directly affect these composite properties are accounted 
for via assumed probabilistic distributions. Collectively, the 
results show that the present technique provides valuable 
information about the composite properties and sensitivity 
factors, which is useful to design or test engineers. Further- 
more, the present methodology is computationally more effi- 
cient than a standard Monte-Carlo simulation technique; and 
the agreement between the two solutions is excellent, as shown 
via select examples. 


Introduction 

Ceramic matrix composites (CMC’s) are potential candidate 
materials for the High Speed Civil Transport (HSCT) structural 
components. As such, they will be required to have an assured 
life of several thousand hours. Estimating the reliability of 
these components is quite a complex process and requires a 
knowledge of uncertainties that occur on various scales. Unlike 
conventional materials, CMC’s have properties that display 
considerable scatter because of the uncertainties encountered at 
every level. For example, there are uncertainties or scatter 
inherent in the constituent (fiber, matrix, and interphase) prop- 
erties and uncertainties in the fabrication process that result in 
scatter in the fiber- volume ratio, the interphase thickness and 
properties, the matrix void-volume ratio, and the geometrical 
parameters of the laminate — such as ply thickness and the 
orientation of the ply. All of these must be formally accounted 
for if CMC behavior is to be predicted with assurance that the 
component will have the required reliability during its operat- 
ing life. Consequently, there is a need for analytical tools that 
quantify the uncertainty in the “response” variables while 
taking into account the inherent scatter in the basic or “primi- 
tive” variables. Primitive variables are the constituent proper- 
ties/parameters that participate at the lowest level (e.g., the 
micromechanics level) in defining a global property. The fiber- 
volume ratio and the individual constituent properties such as 
fiber modulus, matrix thermal expansion coefficient, and fiber 
tensile strength are some of the primitive variables. These are 
assumed to be independent and have their own statistical 
distributions. Response variables are those that characterize the 
composite behavior, such as the longitudinal modulus, the 
transverse thermal expansion coefficient, and the longitudinal 
thermal conductivity. 
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In the current practice of deterministic approaches, uncer- 
tainties are usually accounted for by using safety factors. This 
approach can often yield overly conservative designs, thereby 
reducing the potential of many advanced composite materials. 
Currently, work is underway at NASA Lewis Research Center 
to incorporate the probabilistic distribution of material- 
behavior- and fabrication-related parameters into the 
micromechanics and the macromechanics of CMC’s. The 
primary objective of this work is to develop an efficient 
computational design tool that could account for all the uncer- 
tainties in the constituent properties in a more rigorous manner 
and also provide the overall composite properties and their 
probabilistic distributions. Such information could then be 
used to design structural components to meet the necessary life 
requirements. In addition to providing more rigor to the analy- 
sis than the so-called safety factor approach does, such proce- 
dures would enhance interpretation of the experimentally 
measured CMC properties that are scattered over a wide range. 
Furthermore, the procedure would help identify the dominant 
variables, those that most influence a specific response , thereby 
providing guidelines for quality control during the fabrication 
of these materials. Thus, the methodology could be applied not 
only to the design but also to the development of a better material. 

The approach taken in the present effort was to combine the 
CMC analysis embedded in the computer code CEMCAN 
(CEramic Matrix Composite ANalyzer, refs. 1 and 2) with the 
fast probability integration (FPI) techniques available in 
NESSUS (Numerical Evaluation of Stochastic Structures 
Under Stress, ref. 3). The role of CEMCAN was to provide the 
functional relationships (micromechanics and macromechanics) 
that tie the constituent properties to the equivalent composite 
behavior. The role of FPI was to perform probabilistic analyses 
by utilizing the properties generated by CEMCAN. The results 
are cumulative probability distribution functions (CDF’s) for 
ply and laminate properties of the CMC’s. A CDF is the 
relationship defined by the value of a property (response 
variable) with respect to its cumulative probability of occur- 
rence. The probabilistic sensitivities of the response variables 
to inherent scatter in primitive variables are obtained as 
byproducts of the FPI technique. Select examples of CDF’s of 
composite properties (mechanical, thermal, and strength- 
related properties) and the related sensitivity factors for a 
unidirectional SiC/RBSN [0 8 ] composite laminate are pre- 
sented here for illustration. 

Probabilistic Ceramic Matrix 
Composite Analysis 

As previously mentioned, the present work adopts an inte- 
grated approach, which is a synergistic combination of two 
inhouse-developed methodologies. The first methodology is 
concerned with CMC micromechanics and macromechanics. 


The second one consists of an FPI technique that takes into 
account the uncertainties occurring on various scales in a 
composite and computes the cumulative probabilistic distribu- 
tion of composite global behavior. A schematic of the inte- 
grated approach is shown in figure 1 . A brief description of the 
two methodologies in CEMCAN and FPI are given in the 
following sections. 

Ceramic Matrix Composite Micromechanics 
and Macromechanics 

Over the past two decades several micromechanics-based 
computer codes have been developed inhouse for polymer 
matrix, metal matrix, and CMC analyses (refs. 1,4, and 5). The 
computer code that is currently being used inhouse to analyze 
CMC’s is CEMCAN. CEMCAN utilizes a novel and unique 
fiber substructuring technique in conjunction with conven- 
tional micromechanics that is based on a mechanics -of- 
materials approach. Reference 1 describes the usage of the 
code, and reference 2 describes the various equations embed- 
ded in it, along with the theoretical aspects. Because CEMCAN 
has a unique fiber substructuring technique, it offers several 
advantages over conventional unit-cell-based micromechanics 
theories. It more accurately represents micromechanical inter- 
facial conditions and provides much greater detail about local 
stresses. As shown in parts of figure 1 , CEMCAN ’s methodol- 
ogy consists of an incremental synthesis of the properties, 
starting from the constituents — namely, the fiber, the matrix, 
and the interphase — to form a slice. The slice-level properties 
are obtained by using composite micromechanics equations, 
which are represented by simplified closed-form equations. 
The slice-level properties are, in general, equivalent elastic 
properties such as moduli, Poisson’s ratios, thermal expansion 
coefficients, and heat conductivities. From slice level to single 
lamina level and, subsequently, to laminate level, the code 
repeatedly applies classical laminate theory to obtain 
composite-level properties and responses. Given a specific set 
of loads, the code can progressively decompose, as indicated in 
figure 1 , retracing the steps followed during synthesis so as to 
yield laminate-, ply-, and slice-level responses and constituent 
microstresses. The code can predict ply- and laminate- level 
thermal and mechanical properties as well as a detailed descrip- 
tion of the resulting microstresses caused by an applied load. It 
also accounts for nonlinear effects due to material nonlinearities 
as well as those due to local stress redistribution resulting from 
progressive fracture. By analyzing the response of CMC’s, we 
can account for fabrication-related parameters. A more detailed 
description of this methodology can be found in reference 6. 

Probabilistic Simulation 

There are a number of approaches available for obtaining a 
probabilistic response from a set of independent variables and 
the expressions describing the response behavior. Monte-Carlo 
simulation is one such fairly common approach to obtaining the 
CDF’s of composite properties, given the probability 
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Figure 1 . — Integrated probabilistic ceramic matrix composite mechanics approach. 


distributions of constituent properties, which are considered as 
independent variables. In this technique, randomly selected 
values of the input variables, which are based on their known 
probabilistic distributions, are used to deterministically com- 
pute the value of the response variable. This is repeated, usually 
several hundreds of times, to build the response probabilistic 
characteristics. In essence, this technique requires a large 
number of simulations to generate CDF’s of output variables. 
Although inherently simple, the large number of output sets 
that must be generated to obtain a reasonably accurate CDF of 
output variables becomes its obvious disadvantage. Further- 
more, if the deterministic computation of the response is 
complicated and time-consuming, the computational costs 
could become prohibitive. Obviously, to circumvent this 
computationally very expensive and time-consuming proce- 
dure, more efficient approaches and algorithms are needed. 

For more than a decade NASA Lewis has been engaged in 
developing efficient probabilistic methods. As a result of this 
intensive program, an FPI (ref. 3) was developed to solve a 
large class of engineering problems. In the following para- 
graphs we describe the problem to be solved and how the FPI 
was used to investigate the probabilistic behavior of 
ceramic composites. 

Let us say that there are n random variables in a problem and 
that we want to use probabilistic analysis to compute the 
probability of occurrence of a certain response function 


Z(X) = Z(Xj,X 2 X n ) (1) 

where Z represents the response variable and X represents the 
random variable. Our aim is to compute the probability that Z 
will have a value less than or equal to a given magnitude Z Q . To 
achieve this goal, the performance function, which describes 
how the mechanics of the system behave, can be cast as a limit 
state function g(X), which can be described as 

g(X)=Z(X)-Z 0 (2) 

Here the objective would be to compute P[g(X) < 0]. In the case 
when Z Q describes a limit indicating failure, g(X) is called a 
failure ftinction, or in a classical sense, the probability of the 
stress exceeding the strength. 

Given the probability density function/^*) of the limit state 
function g(x) 9 we can formulate the limit state probability 
P[g < 0] as 

p f = ffeW ^ °] = J • ■ • J* /* M dx (3) 

n 

where Q describes the domain of integration. This multiple 
integration is, in general, very difficult to evaluate analytically. 
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However, FPI has been found to be an excellent tool to evaluate 
equation (3) efficiently and accurately. 

FPI is a probabilistic analysis tool that implements a variety 
of methods for probabilistic engineering analysis and design 
(see schematic in fig. 2). In general, FPI requires the following: 

(1) The independent and uncorrelated input (design) vari- 
ables and their probability distributions must be defined. 
Constituent properties, FVR, void-volume ratio (VVR), ply 
thickness, ply alignment, and the like must be independent 
variables that determine the composite properties. In equation 
(1), X represents these variables. 

(2) There must be a function (called the performance func- 
tion) that defines the relationship between the response variable 
and the independent variables. Ply or laminate properties are 
the response variables in this paper. In equation ( 1 ), variable Z 
is a dependent variable whose uncertainties must be computed. 

The variables that are uncertain in nature are identified as 
independent in step (1). Their probability distributions can be 
obtained from the available measured data or can be assumed 
on the basis of experience and judgment. For most problems, it 
is difficult to determine analytical expressions representing a 
relationship between independent and dependent variables. In 
the case of CMC’s, it is very complicated to build relationships 
for the ply or laminate properties as a function of constituent 
properties, fabrication parameters, and such. CEMCAN code 
has a built-in capability to form such relationships, using 


micromechanics and macromechanics theory. A performance 
function is developed by using a numerical approach. In this 
approach an explicit response function is developed by perturb- 
ing the independent random variables about their magnitude 
and using CEMCAN to compute the response. Discrete evalu- 
ations of the response variables for the perturbed values of 
independent variables are then fit into a function by using 
regression analysis. The uncertainties of a response variable are 
quantified in the form of a CDF by the following procedure: 

( 1 ) The primitive variables and the corresponding probabilis- 
tic distributions are selected. (For example, to generate the CDF 
of the composite longitudinal modulus, the primitive variables 
could be the fiber modulus, matrix modulus, FVR, and so 
forth.) For a given set of values of primitive variables, the 
micromechanics and macromechanics in the CEMCAN com- 
puter code are used to compute the desired response variable. 

(2) The whole process is repeated to generate a table of 
response variable values that correspond to the perturbed 
values of the primitive variables. 

(3) The FPI analysis then uses the previously generated table 
to compute the CDF and the corresponding sensitivities of the 
response. 

FPI computes the CDF of the performance and the proba- 
bilistic sensitivity factor. The CDF can be computed at user- 
defined performance levels or probability levels. Several 
methods can be used to compute a CDF and/or failure probability. 


CEMCAN performance 
function z = f{x-| , x 2 , * 3 ) 


statistics 



Distribution type 


Random variable 


r 



Sensitivity factors 


Fast probability 
integration (FPI) 
analysis engine 



— — _______ 

Output options 




Response cumulative 
distribution function 


(CDF) 


Figure 2. — Fast probability integration input-output schematic. 
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In addition to the CDF’s of the response, the FPI technique 
provides additional information regarding the sensitivity of the 
response with respect to the primitive variables. Sensitivity 
information is very useful for studying the probabilistic varia- 
tion of the response. Note that the sensitivity obtained here 
should not be confused with the conventional deterministic 
sensitivities of the response. Sensitivity in probabilistic analy- 
sis consists of a product of two parts. The first part has the usual 
meaning, akin to deterministic sensitivity. The second part is 
the result of the scatter in a primiti ve variable, and it accentuates 
the deterministic sensitivity. The magnitude of the sensitivity 
factor provides a way to rank the importance of the individual 
physical variables. The sign of the sensitivity factor indicates 
how a specific random variable influences the failure probabil- 
ity. For example, a positive sensitivity value indicates an 
increase in failure probability with an increase in the random 
variable. Thus, the sensitivity information obtained from FPI is 
very useful from the design point of view. The scatter of the 
significant variables can be controlled to improve the reliabil- 
ity. Thus, weak physical variables with large uncertainties may 
have probabilistic sensitivity factors that are more important 
than strong physical variables with small standard deviations. 
Variables with no scatter (deterministic) will obviously result 
in zero values for the sensitivity, which implies that the response 
scatter is unaffected by such variables. 


Results and Discussion 

The material system chosen for this work was SiC/RBSN 
(silicon carbide SCS-6 fiber in a reaction-bonded silicon nitride 
matrix). This composite system is known to have a very compli- 
ant and weak interphase. The mean values of the constituent 
properties are shown in table I, which is taken from reference 7. 
The uncertainties in the composite system occur at various 
levels owing to both the inherent scatter in the material proper- 
ties and the uncertainties associated with the manufacturing 
processes. They may occur at constituent (fiber, matrix, or 
interphase) level, at ply level (FVR, VVR, thickness of the 
interfacial region, etc.) and at composite level (ply angle, layup, 
etc.). 

However, in this preliminary work, uncertainties in the 
constituent material properties and those at ply level were 


considered to quantify uncertainties in some of the composite- 
level properties for a unidirectional [0] 8 composite. Specifi- 
cally, in this work, the FVR, fiber modulus (Zy), matrix modulus 
( E m ), interphase modulus (Zs .), interphase thickness (/■), coef- 
ficients of thermal expansion of the fiber, matrix, and inter- 
phase (cy, a m , and ey), thermal conductivities of the fiber, 
matrix, and interphase {Kj , K m , and K '■), and the room tempera- 
ture tensile strengths of the fiber, matrix, and interphase (5^, 5 m , 
and 5-) are considered random; other parameters were assumed 
to be deterministic. The interphase is a distinct region with a 
finite thickness that arises as a result of fiber coatings or as a 
“reaction-zone” caused by a chemical reaction between the 
fiber and the matrix material as these composites are processed 
at very high temperature. In any case, interphase material is 
assumed to have its own properties, which are different from 
both fiber and matrix properties. For illustrative purposes the 
assumed distribution types and parameters for the selected 
random variables are shown in table II. However, the actual 
distributions for primitive variables were based on experimen- 
tally measured values whenever available. The micromechanics 
and macromechanics embedded in the CEMCAN computer 
code were used to compute composite properties. The micro- 
mechanics equations are not shown here for the sake of brevity , 
but are readily available in reference 1 . The FPI technique was 
used to obtain CDF’s of the required properties. FPI offered 
very valuable additional information in the form of sensitivity 
factors, which represent the sensitivity of output variable 
uncertainty to the uncertainty in the selected primitive random 
variables. 

Three composite properties, the longitudinal modulus E cW 
inplane transverse modulus ^c22’ and inplane shear modulus 
G c 12’ were chosen to illustrate the probabilistic mechanical 
behavior of the composite. The composite coefficients of 
thermal expansion a cll and a c22 and the composite thermal 
conductivities K c 11 and K c22 were chosen as the composite’s 
representative probabilistic thermal behavior. In addition, the 
probabilistic behavior of unidirectional composite longitudinal 
tensile strength at room temperature was also examined. The 
CDF’s of these response variables, along with their sensitivity 
factors with respect to the appropriate primitive variables, are 
shown in figures 3 to 18 and are discussed in the following 
paragraphs. 

The CDF of composite longitudinal modulus E cU and its 
probabilistic sensitivity to the various primitive variables are 


TABLE I.— CONSTITUENT PROPERTIES OF SiC/RBSN COMPOSITE 4 


Property 

SiC fiber 

RBSN matrix 

Interphase 

Modulus, GPa (Mpsi) 

Poisson’s ratio 

Shear modulus, GPa (Mpsi) 

Coefficient of thermal expansion, K _l (°F _I ) 
Thermal conductivity, W/m-K (Btu/ft-hr °F) 

390 (56.6) 
0.17 
117(17) 

4.1 xlO -6 (2.3x KT 6 ) 
22x10^(12.7x10^) 

110(15.95) 
0.22 
45 (6.5) 

2.9X10 -6 (1.6x 10 -6 ) 
5x 10 -6 (2.9x10^) 

3.5 (0.5) 
0.22 
1.4 (0.2) 

2.0XHT 6 (l.lxlO -6 ) 
2.0x10^ (1.2X10" 6 ) 


d Fiber and matrix properties are taken from reference 7 and interphase properties are based on calibration explained in 
reference 7. 


NASA TM-4766 


5 




TABLE II.— PRIMITIVE INPUT VARIABLES DISTRIBUTION PARAMETERS 


Variable 

Units 

Distribution 

Parameter 1 
mean value* 

Parameter 2 
coefficient 
of variation 

Fiber 

Modulus, Ej r 

GPa (Msi) 

Normal 

390 (56.6) 

0.05 

Coefficient of thermal 





expansion, a ^ 

lO^/K (10“VF) 

Normal 

4.4 (2.4) 

0.05 

Thermal conductivity, K j- 

W/m-K (Btu/ft-hr°F) 

Normal 

22 (13) 

0.05 

Tensile strength, Sj 

GPa (ksi) 

Weibull 

2 (285) 

0.05 

Matrix 

Modulus, E m 

GPa (Msi) 

Normal 

110(15.95) 

0.1 

Coefficient of thermal 





expansion, a m 

10^/K (10~VF) 

Normal 

2.1 (1.2) 

0.1 

Thermal conductivity, K m 

W/m-K (Btu/ft-hr°F) 

Normal 

4.2 (2.4) 

0.1 

Tensile strength, S m 

MPa (ksi) 

Weibull 

93 (13.5) 

0.15 

Interphase 

Modulus, E t 

GPa (Msi) 

Normal 

3.5 (0.5) 

0.15 

Coefficient of thermal 





expansion, a, 

lO^/K (10~VF) 

Normal 

2.1 (1.2) 

0.1 

Thermal conductivity, K l 

W/m-K (Btu/ft-hr D F) 

Normal 

2.4 (1.4) 

0.1 

Tensile strength, 5, 

MPa (ksi) 

Weibull 

80(11.6) 

0.15 

Thickness, t, 

pm (mils) 

Normal 

4.2 (0.17) 

02 

Fiber volume ratio, FVR 

— 

Normal 

0.36 

0.1 


“Data obtained from reference 2. 


shown in figures 3 and 4. The computed mean value of E cU for 
this particular composite is 1 83 GPa with a scatter range of 1 55 
to 224 GPa. In other words, if we were to experimentally 
determine the longitudinal modulus of such a composite sys- 
tem, we could expect values anywhere in the range indicated by 
the scatter. Such information obviously gives a test engineer 
useful insight with which to design and plan the test setup and 
the number of tests for a particular material study. 

Sensitivity factors provide information for ranking the 
importance of various input variables in controlling the output 
response and scatter. For example, JE* cl j is most sensitive to the 
FVR (fig. 4). This is to be expected because the longitudinal 
modulus is essentially controlled by fiber-dominated behavior. 
Also important are the matrix modulus, fiber modulus, and the 
interface thickness. If we wish to control the scatter in the 
longitudinal modulus of this particular composite, the biggest 
payoffs will result from controlling the scatter in the FVR. On 
the other hand, the modulus of the interphase material has no 
effect on the scatter of the ply longitudinal modulus. In this 
situation, any effort in the processing or elsewhere that will 
change the modulus of the interphase will not help control the 
scatter in longitudinal modulus. Another point of interest is that 
the sensitivities remain constant throughout the probability 
range considered. In order to verify these results, a limited 
number of Monte-Carlo simulation studies were also under- 
taken. The CDF for the longitudinal modulus obtained from 
10 000 samples is also shown in figure 3. As can be seen from 


the figure, the Monte-Carlo results agree very well with those 
from FPI. Sensitivity information based on the Monte-Carlo 
simulation was not obtained because doing so would have 
involved a substantial additional effort with little added value, 
and it is beyond the scope of the present study. 

Figures 5 to 8 show results pertaining to the inplane trans- 
verse modulus E c 22 and the inplane shear modulus G cl2 . These 
figures also show the range of values that these moduli can take. 
However, the sensitivities do not appear to remain constant 
throughout the range of probabilities. Consequently, impor- 
tance ranking of the primitive variables can change, depending 
on the probability level. For example, in the case of 22 
(fig. 6), the matrix modulus is the most sensitive primitive 
variable at low probability levels, whereas at higher probability 
levels, the interphase thickness becomes the most sensitive 
parameter — a fact not apparent from the micromechanics equa- 
tions alone. The trends of the curves for the transverse modulus 
and the inplane shear modulus appear to be similar in nature. 

Figures 9 through 12 show similar plots for the coefficients 
of thermal expansion a cU and a c 22 - For a cll the predicted 
mean is 3.38, and the scatter range is 2.85 to 3.91 ppm/K. The 
scatter range is based on three standard deviations from the 
mean. An examination of the sensitivity factors indicates that 
the most important variable for controlling response scatter is 
the thermal expansion coefficient of the fiber. Also, the coeffi- 
cient of thermal expansion of the matrix, the FVR, and the 
matrix modulus appear to have a significant effect on the 


6 


NASA TM-1766 



























Cumulative probability 



125 145 165 185 205 225 245 


Composite longitudinal modulus, E c11 , GPa 

Figure 3. — Cumulative distribution function of longi- 
tudinal modulus E c11 of [0] 8 SiC/RBSN composite. 
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Figure 4. — Sensitivity factors for E c11 for [0] 8 SiC/RBSN 
composite. 



30 60 90 120 


Inpiane transverse modulus, E c2 2 , GPa 

Figure 5. — Cumulative distribution function of inplane 
transverse modulus E c2 2 of [0] 8 SiC/RBSN composite. 
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Figure 6. — Sensitivity factors for inplane transverse 
modulus E c22 of [0] 8 SiC/RBSN composite. 



Inplane shear modulus, G c12 , GPa 

Figure 7. — Cumulative distribution function of inplane 
shear modulus G c1 2 of [0] 8 SiC/RBSN composite. 
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Cumulative probability 



Probability 

Figure 8. — Sensitivity factors for inplane shear modulus 
G c1 2 of [0] 6 SiC/RBSN composite. 



Longitudinal thermal expansion coefficient, 
<*c1 1 1 Ppm/K 

Figure 9. — Cumulative distribution function of coef- 
ficient of longitudinal thermal expansion a c11 of [0] 8 
SiC/RBSN composite. 



0.0 0.2 0.4 0.6 0.8 1 .0 

Probability 

Figure 1 0. — Sensitivity factors for coefficient of longi- 
tudinal thermal expansion a c1 i of [0] 8 SiC/RBSN 
composite. 



1.00 1.42 1.83 2.25 2.67 3.08 3.50 

Transverse thermal expansion coefficient, 
a c22> PPnVK 


Figure 1 1 . — Cumulative distribution function of coefficient 
of transverse thermal expansion ot C 22 °f [OJs SiC/RBSN 
composite. 
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igure 1 2. — Sensitivity factors for coefficient of trans- 
verse thermal expansion a C 22 of [0] 8 SiC/RBSN 
composite. 



Figure 13. — Cumulative distribution function of longi- 
tudinal thermal conductivity K^-j of [0] 8 SiC/RBSN 
composite. 


response scatter. The remaining variables have only minimal 
influence on the response. This behavior is to be expected since 
the longitudinal thermal expansion of a unidirectional compos- 
ite is essentially controlled by the fibers. The expressions for 
the thermal expansion coefficients are, in general, much more 
nonlinear because they involve products of moduli, volume 
ratios, and constituent thermal expansion coefficients. We, 
therefore, decided to check the FPI results with a standard 
Monte-Carlo technique as well, and generated 1 0 000 samples 
for the Monte-Carlo simulation. As shown in figure 1 1 for the 
transverse thermal expansion coefficient, the CDF calculated 
by using FPI agrees very well with that determined by the 
Monte-Carlo method. The predicted mean for a c2 2 is 2.25 ppm/ 
K, and the scatter range is 1 .58 to 2.92 ppm/K. The sensitivity 
plot (fig. 12) shows that the dominant variables for a c22 
response are different from those of a cl] . Of the sensitivities, 
the matrix coefficient of thermal expansion influences the 
response probability most, followed by the interphase thick- 
ness, the interphase modulus, the matrix modulus, and the fiber 
coefficient of the thermal expansion. The fiber modulus, which 
was the most important variable for the case of , has little 
influence on the a c22 response. Likewise, the interphase ther- 
mal expansion coefficient has little influence on the response 
scatter. These results are, to some extent, intuitively obvious. 
The behavior in the transverse direction is essentially matrix 
dominated, and therefore, the matrix properties have a greater 
effect on the a c22 response scatter. From these two figures we 
can conclude that the scatter in the transverse thermal expan- 
sion coefficient can be effectively reduced by reducing the 
scatter in the matrix thermal expansion coefficient and the 
interphase thickness. 

The curves for the longitudinal and the transverse thermal 
conductivities are shown in figures 13 through 16. The pre- 
dicted scatter range for K cl \ is 7.4 to 1 1 .6 W/m-K and for K c22 
is 4.2 to 7. 1 W/m-K. The FVR is the most significant variable 
for K cl 1 , whereas for K c22 it is the matrix thermal conductivity. 
All the primitive variables considered appear to significantly 
influence the scatter in the longitudinal thermal conductivity 
response. However, the thermal conductivities of the fiber and 
the interphase appear to have only minimal influence on the 
transverse thermal conductivity response. 

A limited study has also been undertaken to assess the 
probabilistic strength of the composite. The room temperature 
ultimate longitudinal tensile strength was investigated fora [0 8 ] 
SiC/RBSN laminate in which stress is redistributed with pro- 
gressive fracture. The results are shown in figures 17 and 18. 
The CDF of the strength response, as predicted by the FPI, 
shows that we could expect a scatter range of 472 to 850 MPa 
with a mean of 661 MPa. The variable that most influences the 
scatter of the strength is the FVR, which is followed by fiber 
strength. Since a higher scatter was assumed in the FVR 
(10 percent) than in the fiber strength (5 percent), the ranking 
may have been influenced. Matrix strength appears to have the 
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Figure 14. — Sensitivity factors for longitudinal thermal 
conductivity of [0]g SiC/RBSN composite. 



Transverse thermal conductivity, K c 22 . W/m-K 

Figure 15. — Cumulative distribution function of in- 
plane transverse thermal conductivity of [0]q 
SiC/RBSN composite. 



Figure 16. — Sensitivity factors for transverse thermal 
conductivity K c 22 of [0]g SiC/RBSN composite. 



Longitudinal tensile strength, MPa 


Figure 17. — Cumulative distribution function of longi- 
tudinal tensile strength S c11 of [0]g SiC/RBSN 
composite. 
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Figure 18. — Sensitivity factors for longitudinal tensile 
strength S c n of [0] 8 SiC/RBSN composite. 


least influence on strength. This would be expected in a brittle 
system since the matrix starts cracking at low stress levels. With 
an increase in the load, the matrix becomes fully saturated with 
cracks and sheds all its load to the fibers. The ultimate tensile 
strength, therefore, does not depend on matrix strength. Matrix 
strength influences the onset of the nonlinearity in the stress- 
strain behavior though. Thus, in order to control the scatter in 
longitudinal strength, we should exercise stricter quality con- 
trol on the FVR and the fiber strength. Both of these issues are 
related to the fabrication process. 

Summary of Results 

A methodology that provides a formal means to quantify the 
scatter in composite behavior has been described in this paper. 
It takes into account the scatter in constituent properties and 
fabrication-related parameters. This methodology combines 
composite micromechanics and macromechanics and a fast 


probability integration technique to generate cumulative distri- 
bution functions for mechanical, thermal, and strength -related 
properties of a unidirectional SiC/RBSN composite. Collec- 
tively, the results show that the present technique can provide 
quantitative as well as qualitative information that can be used 
as a guide in the fabrication and testing of the material. The 
present methodology is computationally more efficient than 
the standard Monte-Carlo simulation technique, and agreement 
between the two is excellent, as is shown in the results. The 
probabilistic laminate properties and responses can be used to 
design a more rigorous component, thereby avoiding the usual, 
overly conservative, conventional deterministic safety factor 
approaches. The results provide valuable information on com- 
posite properties and sensitivity factors, which is useful to 
design, test, and manufacturing engineers. 
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